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Abstract. We show how the Dixon's system of first order equations of motion 
S . ' for the particle with inner dipole structure together with the side Mathisson 

constraint follows from rather general construction of the 'Hamilton system' 
developed by Weyssenhoff, Rund and Crasser to describe the phase space coun- 
terpart of the evolution under the ordinary Euler-Poisson differential equation 
of the parameter-invariant variational problem with second derivatives. One 
concrete expression of the 'Hamilton function' leads to the General Relativis- 
tic form of the fourth order equation of motion known to describe the quasi- 
■ classical 'quiver' particle in Special Relativity. The corresponding Lagrange 

function including velocity and acceleration coincides in the flat space of Spe- 
cial Relativity with the one considered by Bopp in an attempt to give an 
approximate variational formulation of the motion of self-radiating electron, 
when expressed in terms of geometric quantities. 
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1. Introduction 

Consider a quite popular and fairly general Dixon [T] system of first order ordi- 
nary differential equation^ 

P '° = -2 R ^ PV± " Spv S afi + S 0a = O, (1) 

S' a fj = P a X fj - PpXa 

written in terms of the covariant derivatives, denoted from here on by prime. 

In the theory of General Relativity such equations should hold along the world 
line of a quasi-classical particle endowed with the inner angular momentum (said 
'spin') S a/3 , responsible for its dipole structure. 
| Among several additional side conditions needed to make system ([T]) solvable 

(see [2]), we choose to focus on the one preferred by Mathisson [3] 

i p S pa = 0. (2) 

Now imagine that someone wishes to construct a sort of 'Hamilton' picture of 
the system (JXJ) under the imposed constraint ((2]). There exists a non-conventional 
approach to do this along the following guidelines. First, try to eliminate the 
variables S a @ by means of taking subsequent differential prolongations of (JU |5J|. 
Further, try to find a variational problem with higher derivatives for thus obtained 
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equations, perhaps, under different constraints. Then pass to the corresponding 
Hanhlton-Ostrohrads'kyj counterpart in terms of the generalized momenta. As the 
last step, compose some geometric quantities S a ^ from the canonical variables: the 
momenta and the velocities. If successful, one regains the system (JXJ) , with the 
constraint ([2]) already satisfied identically. 

We show to the end of this paper that (QP follows from a fairly general setting 
of the second order parameter-invariant variational problem as its 'Hamiltonian' 
counterpart by the appropriate definition of S a @ . 

In flat space-time of Special Relativity the differential elimination of the vari- 
able S a @ from ([I] [2]) leads to the fourth order equation of motion 



m 



a 



k 2 - — )x = 0, (± • x) = 1 , (3) 



where k 2 — (x ■ x) is the first integral of ([3]), and 



° a = M ea0pvUf>SPV ■ (4) 

Equation ((3|) was shown by Riewe [4] and Costantelos [5] to describe 'Zitterbewe- 
gung' (quiver) of a quasi-classical particle. 

We show to the end of this paper that |3p occurs as the natural parameterisation 
of an Euler-Poisson (said variational) equation constrained to the manifold k = ko, 
of some parameter-invariant variational problem of the second order. 

The above programme for the flat space-time was carried out in two preceding 
papers [3 7J0 

In present paper rather that go all way round the procedure mentioned above, 
we merely offer a straightforward generalization oh the 'Hamiltonian' depiction 
obtained in [7] to the case of (pseudo)Riemannian geometry. 

2. The Grasser-Rund-Weyssenhoff canonical equations. 

In the space of the fourth order Ehresmann velocities T 4 M let us stick to the 
commonly recognized coordinates x = {x a } € M, u = x = 4^(0), u = 4-f (0), 
u = j-f(0), u = ^1(0). A function £(x,u,u) defined on T 2 M, constitutes a 
parameter-invariant variational problem 5 J Cdr = if and only if it satisfies the 
now well known Zermelo conditions: 

u a -^- = (5.1) 

dii a y ' 

u a ^+2u a ^-C = i). (5.2) 
du a dii a v ' 

We also recall the definition of the Legendre transformation, that is the mapping 
Le: T 3 M T*(TM) over TM given by 

P « = |, (0.1) 
Fir 

P=£-2W 1} , (6.2) 



2 A technical mistake that slipped in the expression for the Hamilton function in paper [6] has 
been corrected in paper 0. 
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where 

d . d .. d 
T> T = u— + u— + u— 

dx du au 

denotes the operator of total derivative, and the canonical coordinates in T (TM) 
are denoted by x, u, p, p^\ Applying T> T to (|5. Ill immediately gives that in terms 
of the mixed set of variables {u,p,p^} the Zermelo conditions look like 

Z 1 ^«V 1) a=0 (7-1) 

Z 2 A = u a p a +u a pW a -£ = 0. (7.2) 
The standard Liouville form A on T*(TM) reads 

A = p.dx + p^- 1 ' .du . 

The system of the canonical equations developed in the paper of Grasser [8] , who 
took as a basis the works of Rund (§] and Weyssenhoff [10] . follow from the exterior 
differential equation 

Le-Hxdk = -XLe~ x dH - f iLe^ 1 dZ 1 . (8) 

In this equation Le" 1 denotes the inverse image operation, acting on forms, with 
respect to the mapping Le, and arbitrary functions A and /i are defined on T 3 M . 
If restricted to the first one of the Zermelo conditions (|7.1|) along the Legendre 
transformation, 

u.pW=0, (9) 

the exterior differential equation JS} defines the Legendre transformation itself, 
along which the function H keeps being constant. It also produces the Euler- 
Poisson equation of the fourth order that demonstrates the parametric ambivalence 
to any local transformation of the independent variable r. 

Let a vector field X on T*(TM) along some curve (x(r), u(t) , p(r) , p^ (r)) be 
its velocity field, 

_ dx a d du a d d Pa d dp^ a d 
dr dx a dr du a dr dp a dr dp\ l ) a 

Then in the coordinate expression the exterior differential equation (|8|) amounts to 
the following system of the first order differential equations [8]: 

Le (10.1) 
oLe + ^u (10.2) 





dx 


dp 




aW = 




du 

aW = 


an 


dp 


o Le — 


x dn o 

dx 


dr 


dr 


o Le — 


\ m o 
du 



(10.3) 

MP W. (10.4) 



3. Implementing the covariant derivation. 
Let us introduce the following change of local coordinates in T 2 M: 

{x, u, u} — > {x, u, u'} , 



4 



R. YA. MATSYUK 



where 

u'^^+rv^. (ii) 

Inspired by ©J we also may consider the 'covariant momenta' given by 

p = a(£or 1 )_ tp „ ) , (122) 



In (|12.ip we have 

dtCoch- 1 ) dC 
= — r ° 

by virtue of (fTTj) . In (112. 2p again on the strength of (jTTJ) one computes 

arro^- 1 ) a/: , , mo . ac , , 



(13) 



(14) 



On the other hand, the rule for the covariant derivative of a covariant vector says: 

OT 



Relation (fT3")) should be understood in terms of the notation (|6.ip as 

ro'^/dopWof 1 , (16) 
from where it immediately follows that also 

ir = ir°^ • (17) 

Inserting ^TBJ) and dTTJ) into (JT5J) and then together with {HI and ([6~Tj) into (fT2~2]) 
gives 

B= (^-^ (I) »)° fl - r ^° fl 

by the definition (|6.2I) . This suggests the corresponding change of coordinates in 
the manifold T*(TM): 

{x,u,p,p w } A {x,M,7r,7r (1) } , 

\ tt q = Ptt of i -rvV 1) p°r 1 . " Sl 

Thus the Legendre transformation Le is represented in the coordinates {x, u, v! , u"} 
and {x, u, it, 7i"W} by (TT2"|) , which is the local expression for 

Z^ = $oLeo0- 1 . (19) 

Proposition 1. Let $j — H o depend on x, u, n, ttW through the invariants 

J = U-U, 1p = 7T.U, = 7T • 7T*- 1 - 1 (20) 



HIGHER ORDER VARIATIONAL ORIGIN OF THE DIXON'S SYSTEM 



5 



only. Then the 'Hamilton equations' i!0\) take the shape 
dx 

- = u ^ (21.1) 

u ' = 2mW(^) ^-oLe + flu (21.2) 
\oip J or) 

7r' a oZe= -R a f} p v u p vPm( 1) v (21.3) 

^'Le = -2u (f^Y^Te-w-faW, (21.4) 
\oip J cry 

where fi and /i from H0\) are related by fj, = [i o . 

Proof. First we compute the derivatives of 7, -0, 77. As far as -0 is a mere contraction 
of the covariant vector n with the contravarient vector u, it contains no metric tensor 
g a p; thus 

. . , (22) 

From the Riemannian geometry we recall the formulas for the partial derivatives of 
the metric tensor 

-q-^T = ffapl ' up + gppt ' va , 



dg af) 

so that 



dx" ' 9 pv 9 f 



Or) 

From flU and (fill) we have 



(23) 



=7^ 24 ' 2 

Moreover, recalling our assumption about the dependencies of .f) upon its argu- 
ments, from (|20[) we get 

=2^)", (25.1) 



dirt 1 ^ dr/ 

ft?=V Utt+ ^ a ' (25 - 2) 



and further on, by the use of (|2"2"]l and 



^ = 2^ tt ,-2|VW (25.3) 
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Again from (|18l) . and recalling the rules for the covariant derivatives of vec- 
tors (jTTJ) and covectors, 

7r' a = ^ - T" Pa ny , (26) 

we obtain 

- r^r^t/W 1 ^ + r^i^V^ 1 ),, . (27.2) 

In view of (|24.2p equation (| 1 . 1(1 now becomes 

^ = A ^2 $ oie . ( 28) 
dr on 

As far as 

^--^u (29) 
and in view of (|28|) . equation (|10. 1[) transforms into 

— = Au— o $ o Le, 
dr dip 



so the choice 

dip) 

seems legal. Thus in the sequel we implement the definition 

dx 



o $ o Le (30) 



(31) 



Let us consider equation (I10.2[) . In it we substitute (|24.1[) together with (|25.ip 

for — 7TT- and afterwards we put A—— o $ o Le = u on the strength of (|2"51) and 

op*- 1 ' on 

of (15T)) . But in view of (|3"U)) this exactly produces (|21.2p . 

Now let us turn to equation (|10.4p . We apply formula like (|2"6")l into (|27.1I) in 
order to use it in (|10.4p together with (|24.3[) and (|25.2p . And one more time we 
use (|2"8")) accompanied by (j3"Tj) and then apply (j3"0")) . This amounts to (|21.4p . 

The equation (| 10 . 3[) is the most interesting. This is the evolution equation. 
Recall that the left hand side there is given by (|27.2I) . in where ^ should be 
substituted by (|26p as usual. 

As the next step for n^' o Le in (|10.3|) we substitute (j21.4j) and for u' in there 
we substitute (12T2"]) . 

In the right hand side we implement formula (124.41) with subsequent use of (|25.3I) , 

and, moreover, inserting there — — from (129p together with (1301) . 

On 
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The final step consists in grouping the remaining terms to fit in the well known 
definition of the Riemannian tensor, which in contraction with u and reads 

n v a n u u ~ n v a r» u u 

oxP ox 

= ^ l KR a p P v Wu^ (32) 

This completes the proof. 

Proposition 2. Equations jTJ 0) follow from \21.!J\) if one introduces 

S a p = u a ir {1) p - upTT {1) a . (33) 
Proof. If S = u A then, taking n^' from (f2"Tj) . one computes 

S' = U A TT^ 1 ' + U A TV = /ilt A 7T^ -llA7T-/lIlA 7T^ = 7T A It . 

Thus we reasonably identify P with tt and this way regain the second equation 
of ©. 

Equation ([2]) is satisfied automatically. 

On the basis of skew-symmetric properties of the curvature tensor one rewrites 
the equation (|21.3p as follows 

n' a = P' a = -Rap^UpTT^u 

= -\V VufiS P v by©, 

and thus regains the first equation of ([I}. 

Remark. Up to this point any concretization of the expression for the 'Hamilton 
function' f) does not matter. 

4. The 'Zitterbewegung' and electron self radiation in covariant 

Riemannian framework 

Let at last the 'Hamilton function' take the following concrete shape: 

^ = 7r . u + M! 7r (i). 7r (i)_A|| u |i. 

4 



Equations (j2l~2"]l and (|2TlS|) look like 



u'= ]p 3 ' 2 uj^+fLu (34) 

Q A 

t/ 1 )' o Le = 7 1/2 (ri o he) u - w - ft + u. (35) 

4 v/7 



Now let us apply to ([34)1 . along the Legendre transformation, the Zermelo condi- 
tion ([9]) in terms of the covariant variable m^ 1 ', given by ([T6]) . This defines jl: 
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To solve the system of equations (f34| l35j) for w, it suffice to differentiate (f34|) and 
to substitute there (l3"5t for w^' and ([M]) for w^ 1 ' again. This resulting in 



— 



6- 



. (u' • uf 



A 



(37) 



As far as we remember that the equation (|21.3I) keeps the property of the am- 
bivalence to arbitrary transformations of the evolution parameter r, we absolutely 
may pass to the natural parameter s for which u ■ u = 1 together with v! ■ u = and 
u" -u = —u'-vf. Then (|21.3p with ir given by (j3"Tj) finally amounts to the following 
dynamical equation of motion: 



D 

ds 



D 2 x D 2 x 



ds 2 ds 2 



A 



Dx a 
ds 



D 3 x a 
ds 3 



(38) 



Proposition 3. In flat space-time equation (38\) on the constraint manifold k = ko 
reduces to the Riewe-Costantelos equation f3|) of the quasi-classical 'Zitterbewegung' 

with the frequency y f^o — 4 ^ v P u ^ n 9 A = fc§ + 2^j- with a from Q). 



given 



Proof. Take respect for the definition k = • Q-r = ko 2 in 

The Lagrange function for the forth order equation (|2f .3|) with w and 
by (32), HMD, and (011), is: 

£ fc = (fc 2 +A)|M| . (39) 

In 1946 Bopp [ 1 1| in the attempt to give an approximate variational formulation 
of the self radiating electron in flat space-time of Special Relativity, introduced a 
Lagrange function that might have been even at that time expressed in terms of 
the geometric quantity k — the first Frenet curvature of the radiating particle's 
world line. The Bopp's Lagrangian is nothing but (f39|) . seen in pseudo-Euclidian 
coordinates. 
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